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Abstract 

Two approaches to describe the thermodynamics of a subsystem that interacts with a thermal 
bath are considered. Within the first approach, the mean system energy Eg is identified with 
the expectation value of the system Hamiltonian, which is evaluated with respect to the over- 
all (system+bath) equilibrium distribution. Within the second approach, the system partition 
function Zg is considered as the fundamental quantity, which is postulated to be the ratio of the 
overall (system+bath) and the bath partition functions, and the standard thermodynamic rela- 
tion Eg = —d(lnZs)/d(3 is used to obtain the mean system energy. Employing both classical 
and quantum mechanical treatments, the advantages and shortcomings of the two approaches are 
analyzed in detail for various different systems. It is shown that already within classical mechanics 
both approaches predict significantly different results for thermodynamic quantities provided the 
system-bath interaction is not bilinear or the system of interest consists of more than a single 
particle. Based on the results, it is concluded that the first approach is superior. 
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I. INTRODUCTION 



Let us consider a system in equilibrium at fixed volume and temperature described by a 
canonical ensemble. If the partition function Z for the corresponding canonical distribution 
is known, all thermodynamic quantities, such as the internal energy, the entropy, and the 
specific heat can be calculated by simple differentiations of In Z with respect to the temper- 
ature T. In many situations it is of interest to consider a subsystem of the overall system, 
e.g. a smaller system that interacts with its environment (in the following referred to as the 
bath). The question arises, how to describe the thermodynamics of the subsystem? 

If the system-bath coupling is weak, then (to leading (zero) order in the system-bath 
coupling) the system can be described by a canonical distribution determined by the corre- 
sponding system Hamiltonian. This finding builds the basis of standard statistical thermo- 
dynamics. If the coupling is not weak, however, the interaction with the bath affects the 
system density matrix, the system partition function Z$, and all thermodynamic quantities 
2I i i 1 1, Furthermore, new interesting effects arise due to the system-bath entanglement 
5, 6, 7,y, 9), Q. In such a case, it is not straightforward, e.g., to define the energy of 
the system unambiguously. Different definitions are possible which incorporate, to a certain 
extent, the system-bath interaction into the system energy. 

There exist two major approaches to describe the thermodynamics of a subsystem 
(strongly) coupled to a bath. The first approach (in the following referred to as approach 
I) considers the mean energy of the system as fundamental quantity and assumes that it is 
given by the expectation value of the system Hamiltonian which is evaluated with respect 
to the total (system+bath) canonical equilibrium distribution The second ap- 

proach (approach II), on the other hand, is based on the partition function of the system, 
Zg, and postulates it to be given as the ratio of the total (system+bath) and the bath par- 
tition functions. Based on the thus defined partition function of the system, the standard 
relations of thermodynamics are invoked to calculate the internal energy, the entropy, and 

nnn 

the specific heat by differentiations of lnZ^ with respect to the temperature T |14 1151. 1161] . 



In recent work 
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131 ]. it was shown that the two approaches give different results for 



the specific heat of a quantum mechanical point particle or a harmonic oscillator bilinearly 
coupled to a harmonic bath, despite the fact that the results are identical in the classical 
case. In the present work, we analyze the two approaches and assess their validity in detail 
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for several different systems including several quantum and classical point particles and 
nonlinear system-bath coupling. 

Our main findings are summarized as follows: If approach I is employed, then the knowl- 
edge of the system partition function Z$ alone is not sufficient to describe the thermody- 
namics of the system and the standard procedure calculating the thermodynamic quantities 
by differentiations of In Zs with respect to T is no valid. Additional knowledge of a bath- 
induced interaction operator As (see below) is required. We derive general expressions for 
the mean system energy, entropy, and the specific heat in terms of Z s and A s . 

If, on the other hand, the strategy of approach II is followed to define the partition 
function of the system Zs, then obtaining the internal energy, the entropy, and the specific 
heat by differentiations of In Zs with respect to T cannot be justified by referring to the 
standard thermodynamic machinery. The so-obtained expressions must be regarded as the 
definitions of the corresponding quantities, whose validity must be proven a posteriori. 

We also show that the differences between the two approaches are not of purely quan- 
tum origin. The approaches predict significantly different results already within classical 
mechanics, provided the system-bath interaction is not bilinear and/or the system of inter- 
est consists of more than a single particle. The general results are illustrated by explicit 
calculations of thermodynamic quantities for several classical and quantum model systems. 
Based on the thus obtained results, we conclude approach I is clearly superior to approach 
II, both on physical and logical grounds. 



II. THERMODYNAMICS OF A CANONICAL ENSEMBLE 

Let us consider a canonical ensemble which is kept in constant volume L at fixed tem- 
perature T. The corresponding equilibrium distribution (density matrix) reads 

p = Z~ 1 exp{-PH}, Z = Tr(exp{-(3H}). (1) 

Here H is the Hamiltonian, Z is the partition function, (3 = l/(k B T) (k B being the Boltz- 
mann constant), and Tr(...) denotes the integration over phase space variables in case of 
a classical ensemble or taking the trace in case of a quantum ensemble, respectively. We 
further define the (information) ensemble entropy operator as 

S=-k B \np = ^{H-F) (2) 
3 



and the free energy 

F = ~lnZ, (3) 
so that the canonical distribution (TjQ) can alternatively be written as 

p = ex P {-(3(H-F)}. (4) 

Averaging Eq. (j2j) over the canonical distribution (j4j) we obtain the following expression for 
the ensemble- averaged entropy 

(S) = ^((H) - F), (5) 
where we have used the notation (...) = Tr(p...). The specific heat is determined as 

Using Eqs. (JI]) and (jlj) we can express the internal energy, the entropy, and the specific heat 
through the free energy: 

W - ~. (8) 

We emphasize that Eqs. (|7|)-([9|) are a direct consequence of the explicit form and 
temperature-dependence of the canonical operator ([TJ. 

The second law of thermodynamics in differential form then reads 

Td (S) = d (H) . (10) 



III. THERMODYNAMICS OF A REDUCED CANONICAL ENSEMBLE 



A. General expressions 

The expressions for the thermodynamical observables of a canonical ensemble listed above 
can be found in any textbook on statistical thermodynamics. To study the thermodynamics 
of a subensemble let us now consider a system that is interacting with its environment (in 
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the following referred to as the bath). We assume that the Hamiltonian H can be split into 
the system (S) Hamiltonian, the bath (B) Hamiltonian, and their coupling, 

H = H S + H B + H SB . (11) 

Here Hs depends solely on the system degrees of freedom and H B depends exclusively on 
the bath degrees of freedom. 

To calculate an observable which depends on the system degrees of freedom only, it is 
expedient to introduce the reduced density matrix of the system, p$, which is defined by 
averaging the total density matrix p over the bath degrees of freedom 

p s = Z~ l Tz B (exp{-f3(H s + H B + H SB )}) . (12) 

Following Kirkwood 171 ]. Eq. (j!2p can equivalently be rewritten as 11] 



Ps = Z s - 1 ex P {-P(H s + A s )}, Z S = Z. (13) 

Here, we have defined the bath-induced interaction operator 

A 5 = -- lnTr B (exp{-(3(H s + H B + H SB )}) - H s , (14) 

which, in general, depends on the system degrees of freedom and on the temperature, T. 

If the system-bath coupling is week (H-ffssll "C \\Hs\\ , then, to the leading (zero) 

order in the system-bath coupling, the interaction operator is given by As = —ln(Z B )/f3. 
Here Z B is the partition function of the bath canonical distribution 

p B = Z^expi-PHs}, Z B = Ti B (ex P {-(3H B }) (15) 

In this weak-coupling case, the overall partition function factorizes, Z = Zs, c Z B (the sub- 
script "c" stands for "canonical" ) and the reduced density matrix (|T2|) is determined by the 
canonical distribution for the system alone 

p s , c = Zsl exp(-(3H s ), Z s , c = Tr 5 (exp{-/3# s }). (16) 

In general, however, Hs B cannot be neglected in Eq. f)14p and the simple expression ffTB"]) is 
not valid. 

Analogously to Eq. (j2J), we introduce the system free energy 

F s =-hnZ s , (17) 
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so that the distribution (fi"3l) can equivalently be rewritten as 



p s = exp{-P(H s + A s - F s )}. (18) 
The distributions (|T3|) and ([18]) can be used for calculating the expectation value of any 



operator Y$ which depends only on the system degrees of freedom 
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(Y S ) = (Ys) s , (-)s = Ttgipg...) = T*s+b(p-). (19) 



B. Approach I (Mean energy approach) 

In approach I, the mean system energy E$ is associated with the expectation value of the 



system Hamiltonian (Hs) [U 
heat is determined as 



3,Q. 



Il3j |. Accordingly, the system contribution to the specific 

To obtain an expression for the entropy of the system, we follow the presentation in Sec. 
2 and define the information system entropy as 

S s = -k B Inps = ^(H s + A s - F s ). (21) 
Averaging the expression over distribution f[T5j) , we obtain the analogue of Eq. §5§ 

(Ss) = ±{{H S ) + (A s ) - F S ). (22) 

In passing, we note that there exists a controversy in the literature whether the so defined 
(S s) c an indeed be considered as a proper thermodynamic entropy of the reduced system 
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191 ] . For the purpose of this paper, Eqs. (12T1) and ( |22l) can merely be considered as 
intermediate mathematical expressions, which allow us to conveniently derive the formulas 
discussed below. 

In contrast to the situation discussed in Sec. 2, the system free energy Fs is no longer a 
universal object which determines all relevant thermodynamic quantities (Ss), (Hs), and Cs 
through Eqs. (J7])-©. It is the explicit dependence of As on the system degrees of freedom 
and the temperature which violates the standard thermodynamic expressions for the reduced 
system. It is straightforward, however, to generalize Eqs. ([7j)-([9]) if we assume that operators 
Hs and A s commute. In this case, using Eqs. (|T7|) - (!20|) . we obtain the expressions 

d(\nZ s ) dF s . , /A , /dA s \ , s 



/dA s \ , d(A s ) , . 

/ \ (26) 



It is important to note that 

\~dT~ / r dT 

due to the explicit temperature-dependence of ps- If H$b ^ 0, then the last two terms in 
Eqs. ( 1251) and (|2"5|) give a non-negligible contribution. If the bath-induced potential A s 
in Eq. ( 12^1) is temperature-independent, then the usual formula (jHJ) for the entropy holds, 
but the expressions for (Hg) fl23l) and (|25p do contain the bath-induced contributions. 
Furthermore, Eq. (|25|) shows that the second law of thermodynamics (in differential form) 
is modified to 

Td (S s ) = d (H s ) - (dA s ) + d (A s ) . (27) 

To summarize, knowing the partition function Z$ or the free energy F$ alone is not 
enough to calculate (Ss), (Hs), and Cs within approach I. Instead, the more general Eqs. 
(I2"3l) -(12T1) must be used. This should be taken into account if, e.g., work theorems are 
employed to obtain the system partition functions beyond the weak system-bath coupling 



limit 
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221 ] . Th e ge neral expressions (I23l)- (|2T1) are also important for the thermodynamics 



of small systems 23|, |2J, |25 |. 



Eqs. ( 123]) - (T27]) have been derived assuming that H$ and As commute. This requirement 
is not as restrictive as it might seem at first glance and is enough for the purposes of the 
present article. It is obviously fulfilled within classical mechanics. It is also fulfilled in 
the semiclassical limit considered below described within the Wigner function formalism, 
provided we start from the total (system+bath) Wigner distribution and introduce the 
reduced Wigner distribution of the system by averaging the total Wigner distribution over 
the bath degrees of freedom (see Sec. IV C). 

Even in the general case, where [Hs, As] ^ 0, analogues of Eqs. (Ej)-(l9]) can be derived. 
To this end, we first rewrite the system density matrix in the form 

p s = Zf exp{-(3H s } exp{-(3A s } (28) 

with a slightly redefined bath-induced interaction operator 

A 5 = ~ In (ex P {(3H s }TiB (exp{-(3(H s + H B + H SB )})) . (29) 



Second, we replace Hs — > 7^5 in Eq. (1281) (with 7 being a numerical parameter) but keep 
As unchanged, so that the distribution fl28|) becomes 

p^s = Z~ s exp{-/?7^}exp{-/?A 5 }, (30) 

Z 7iS = Tr 5 (exp{-/?7# 5 } exp{-/3A 5 }). (31) 

Then Z 7i s and F 7i s = —ln(Z 7j s)/P become 7-dependent and, we obtain the expressions 

(32) 



1 cf(ln Z yS ) dF ljS 



(3 c?7 c?7 



7=1 

(33) 



7=1 



d-fdT 

C. Approach II (Partition function approach) 

The reduced distribution (TT3T) remains unchanged if we introduce a certain (possibly 
temperature-dependent) function T(T) and redefine the system partition function and the 
bath-induced operators as 

Z S ^Z S /T(T), A s ^ A 5 + ilnT(T). (34) 

The analogues transformation for the distribution ffl8l) reads 

F s ^F s -hnT(T), A 5 ^ A s + ilnT(T). (35) 

The transformations (13"4"1) and flHB]) shift the origin of the bath-induced interaction Ag. For 
example, we can take T(T) = (the bath partition function is defined via Eq. ffl5l) ). This 
choice is especially reasonable for the weak system-bath coupling, because it makes A5 = 
if Hsb = 0. In general, for any Hsb 7^ 0, it yields the system partition function 

Z s = Z/Z B . (36) 

The definition (136]) is the key to the approach II. According to the recipe developed in Refs. 
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m 



lq . we should identify fl36l) with the system partition function, and use the standard 
Eqs. (JT])-© to calculate the necessary thermodynamic quantities. 

The approaches I and II result in different expressions for thermodynamic quantities, 
provided the system-bath coupling is not weak. This will be illustrated in Sec. 4 based on 
different examples. The validity and predictions of the two approaches are discussed in Sec. 
5. 



S 



IV. ANALYSIS OF THE TWO DIFFERENT APPROACHES FOR ILLUSTRA- 
TIVE EXAMPLES 



A. Model system 

We consider a general system-bath problem. We assume that the system (e.g., a 
(macro)molecule) consists of Ns point particles, Pi, and Mj being their positions, mo- 
menta, and masses. The bath comprises Nb point particles with positions Xj, momenta p^, 
and masses mj. All interactions (S — S, S — B, and B — B) are pairwise, so that the parts 
of the overall Hamiltonian H = H$ + Hb + Hsb are explicitly written as follows 

N S , p2 >> N S 

Hs = ^2\2M- + Us{Xl) + ^ Uss{Xt ~ (37a) 

i=l ^ 1 ' i>j 

N B 2 N B 

i=l 1 i>j 
N s N B 

Hsb = Y,Y, U SB{Xi - x ). (37c) 

i=l j=l 

Here Uss, Ubb, and Usb are the corresponding interaction potentials (which may be different 
for any pair of particles i and j) and the system is allowed to be subjected to an external 
potential Us- For clarity, we consider a one-dimensional ensemble. A generalization to the 
three-dimensional case is straightforward. 



B. Classical mechanics 



We first consider the two approaches to subensemble thermodynamics for a system of 
classical point particles. In classical mechanics, the partition function for any canonical 
distribution is a product of the momentum and coordinate contributions. Furthermore, the 
momentum contributions to the reduced distribution fflBl can be integrated out, so that the 
only nontrivial part of the distribution is the contribution of the potential energy. 



1. A single Brownian particle 

Let us first consider a single Brownian particle, which corresponds to Ns = 1 in Eq. 
(I37al) . After the insertion of Hamiltonians (I37al) . (I37bl) . and fl37cj) into Eq. (Tl2]) . we can 
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make use of the isotropy of space and change the integration variables Xi — > X{ — X. This 
way, we obtain As = and arrive at the standard result that the reduced distribution (fl3"]) 
is the canonical distribution (TIB]) determined by the system Hamiltonian 

p s = Z s l eM-P(^ + U S (X)}, (38a) 

Z s , c = ^2nM/(3 J dXexp{-(3U s (X)}. (38b) 

It is tempting to assume that all thermodynamic characteristics of the Brownian particle 
can be obtained through the bath- independent, canonical system distribution (j38aj) or, what 
is equivalent, through the differentiation of the partition function (I38bl) according to the 
standard Eqs. ©-([9]). This is indeed the case if we use approach I. If we follow approach, 
however, we obtain 

Z/Z B = V Zs,c. (39) 
Here Zs c is the free-particle partition function (I38bl) and 



77 = ^ f / ^"Tl — ■ ( 40 ) 



J dx\...dxN B exp j 






'} 


J dx\...dx 


7v B exp j 


[-/? (ES ^b^-x,)) 







Apparently, r/ ^ 1, in general. This becomes evident, e.g., if we expand the numerator in 
Eq. (14"0|) in powers oiUgB- Symbolically, t] = 1 + O {\\U sb\\) ■ In approach II, the factor of r\ 
induces (unphysical) bath-dependence of the system mean energy, entropy, and the specific 
heat. Thus, approaches I and II lead, in general, to different predictions even for a single 
classical Brownian particle. 

An important exception is a Brownian particle bilinearly coupled to a harmonic bath, 

/ 2 2 2 \ 



8=1 



N B 2 

ff*B = £^ (X 2 - 2Xx.) , (41b) 

i=l 

where u>i denotes the frequencies of the bath oscillators. In this case, we obtain 77 = 1 and 
the two approaches give the same result. This is only the case for the simple form of the 
bilinear system-bath coupling. 
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If we retain the harmonic bath (14 lap but add a nonlinear interaction term to Hsb, the 
situation differs. Let us consider, for example, the potential 

ff _g{^ (Ml0 + _jy (42) 

where £j denote the corresponding constants. The additional term in the potential, the 
form of which has been chosen for demonstrative purposes, may describe repulsion of the 
particles at short distances. The reduced system partition function Z$ is given by Eq. ( I38bl) . 
Incorporating Eqs. (141 al) and (f4"2"|) into Eq. (I4(jp . we obtain for the factor rj, which describes 
the deviation of the system partition function Z$ from the ratio Z/Zb, 

N B 

V = Hexp{-f3u J ^2m 3 t j }. (43) 
i=i 

Thus, even for this rather simple example, the factor rj can significantly differ from unity and 
also acquire a temperature dependence. Within approach II, this would result in incorrect 
predictions for (Ss), (Hs), and Cs [26j |. 



2. A harmonic dumbbell 



Let us suppose that the system (e.g., a molecule) consists of a collection of point particles. 



In general, the explicit evaluation of the bat 



bath coupling limit is a difficult task H 



i-induced potential As beyond the weak system- 



27 



281 ]. except in the case of a harmonic bath 



(modelling, e.g., a Gaussian solvent) bilinearly coupled to the system. In the latter case, 
the integrations over Xj in Eq. ffl2l can easily be performed analytically. It is instructive to 
consider the simplest nontrivial situation, when the system consists of two identical particles 
(Ns = 2). Such a model can describe, for example, a diatomic molecule or a dumbbell. If we 
require the total (harmonic) Hamiltonian H (fTTl) to be translationally invariant, we arrive 
at the expression 

*8 = 4h + ih+ ^(Xi - X 2 )\ (44a) 



H. 



SB 



N B 

E 

i=l 



2M 

Pi 
2m,- 



+ w4 (Xi 



+ 



nijuo 



{X 2 - Xi)'< 



(44b) 



where uu, cj 2 i, and Qs denote the corresponding oscillator frequencies (29]. 
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For later use it is convenient to rewrite Eq. (144b[) in the equivalent form 
H 3 B + H B -^\— + -(u li + u 2i )^x i - + | \+—( Xl -X 2 ). (45) 



i=l 

Here 



Nr 2 2 

rrijUjt-ujo 



nl = T T P^ T , (46) 

is the frequency of the bath-induced harmonic potential. Integrating out the bath modes, 
we obtain a reduced distribution p s of the form of Eq. f[T3"l) with 

A s = — A( Xl -X 2 ) 2 , (47a) 



_ 2tvM / 2tt 

Zs "TlP(o[Tn|)' (47b) 

Here L is the (one-dimensional) system volume. It is seen that the influence of the bath 
manifests itself in the additional attractive harmonic potential As, which is coordinate- 



dependent but temperature-independent, so that dAs/dT = 32]. 

For the Hamiltonian (144a[) - (144b[) . the ratio of the total and the bath partition functions 
yields the system partition function, i.e. Z/Z B = Zs (Eq. (|47bl) ). One might thus expect 
that approaches I and II give the same predictions for the thermodynamics quantities. Due 
to the presence of the bath-induced potential As, this is, however, not the case. Indeed, if 
we would follow approach II, the mean system energy is given by 

dQnZs) _ d(\n{Z/Z B }) _ 3 

s d/3 dp 2(3' 1 ' 

which corresponds to the thermal energy of a system with three degrees of freedom (one for 
the center of mass translation and two for the vibration). If, on the other hand, according 
to approach I, we associate E s with (H s ), then Eq. (123]) yields 

Thus, approach II predicts for the heat capacity Cs = 3k B /2, irrespective of the strength 
of the dumbbell (0|) and solvent-induced (£1%) potentials. This seems to be physically 
incorrect, given that the reduced distribution ps contains the bath- induced attractive har- 
monic potential As (147aj) and does not coincide with the canonical distribution ( TL6l) for the 
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dumbbell alone. Approach I, on the other hand, predicts the bath-dependent specific heat 
to be 

The coupling to the solvent has the strongest influence in case of two free Brownian particles 
= 0). In this case, the actual number of degrees of freedom is reduced by one, which is 
in accordance with the physical expectations. 



C. Quantum mechanics 



We next analyze the two approaches to subensemble thermodynamics for quantum me- 
chanical point particles. In the quantum mechanical case, there exist pitfalls and subtleties 
in the calculation of the specific heat already for a single Brownian particle bilinearly coupled 



to a heat bath of harmonic oscillators. This has been demonstrated in Refs. 
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13. To 



elucidate the nature of these subtleties and to simplify the presentation, we restrict ourselves 
to a semiclassical analysis and calculate the leading order (~ h 2 ) quantum corrections to 
the thermodynamic quantities. To this end, we employ the Wigner representation 33]. 

To simplify the notation, we introduce a collective index a which runs over all iVg + Nb 
system and bath particles and use a tilde to denote the corresponding positions, momenta, 
and masses. Thus, the total Hamiltonian (1111) reads 

Within the Wigner representation, we treat p a and x a as (semi) classical phase space vari- 
ables. The Hamiltonian retains its classical form, but the canonical distribution for the 
overall system flTJ is given by the corresponding Wigner distribution (denoted by the super- 
script W) 

p w = Z~ l exp{-(3H + A# (1) } + 0(A 2 ), (52a) 

Z = Tr(exp{-/3# + XH W }) + 0(A 2 ), (52b) 

A = (2nh) 2 . (52c) 

Here, as in the classical case, Tr(...) denotes the integration over the corresponding phase 
space variables. The quantum correction is explicitly given as 33|, [34] 



(3 2 d 2 U (3 



+ 



8m a dx\ 24m a \dx a 



dU 




P PaPb d_U_ 
24m a m b dx a dxb 



(53) 
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It is noted that is explicitly temperature-dependent and contains mixed coordinate- 
momenta terms. 

The reduced Wigner distribution of the system, p<T, is obtained by integrating Wigner 
distribution of the overall system, Eq. (I52al) . over the phase space variables of the bath. Since 
H and in Eq. (152aj) are functions but not operators, expressions (jl~2l - (l2Tl) derived in 
Sec. 3 remain also correct for the reduced Wigner distribution. In general, the Wigner 
transform of an operator differs from its respective classical expression, because quantum 
mechanically positions and momenta do not commute. However, if an operator can be split 
into a part which depends only on coordinates and a part which depends only on momenta, 
then the Wigner transform is given by the corresponding classical expression. This is the 
case for H$, Hg, and Hsb- Therefore, we can use the classical expression for H$ while 
evaluating the mean system energy. 



1. A quantum Brownian particle 



As in Ref. Il3j, we consider a single Brownian particle bilinearly coupled to a heat bath 
of harmonic oscillators. The system Hamiltonian is given by Eq. (I37al) with N$ = 1, the 
bath Hamiltonian by Eq. ( 141 aft and the system-bath coupling by Eq. ( 141bl) . Inserting 



the corresponding formulas into Eqs. (|5ip and (153]) we can integrate the bath degrees of 
freedom out of the overall (system+bath) Wigner distribution (152aj) and arrive at the system 
distribution 

P Y = Z~ s l exp{-(3(H s + A S )} + 0(A 2 ). (54) 

Here H$ = P 2 / (2M) is the free particle Hamiltonian, the bath-induced interaction operator 
reads 

3=1 

and the partition function is given by the expression 



2;rU lJ^(1 + A(^) 2 /24)) + 0(A 2 ). (56) 



/5(1 - A(0/5) 2 /12) V P 

The bath-induced operator A5 is position- independent, but depends on momentum and 
temperature. 
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The ratio of the total (Z) and bath (Zb) partition functions can also be readily obtained 
from expressions (I52ap - (l53]) to yield 



Z/Z B = L^^-(l- A(fi/3) 2 /24)) + 0(A 2 ). (57) 

A comparison of Eqs. (1571) and (I56p shows that the partition function of the system is 
not given by the ratio of the total (Z) and bath (Zb) partition functions. The quantum 
corrections in the two expressions have the same magnitude but opposite signs. 

Since the bath-induced operator A s is momentum and temperature dependent, it is not 
expected that the differentiation of the partition function alone gives the averaged energy 
(Hs). Indeed, the calculation gives 

1 AO 2 

- d(\n Z s )/d(3 =^-^P + 0(X 2 ). (58) 

On the other hand, if we follow the approach I and use Eq. (1231) or perform directly an 
average over the distribution ( |54"1) . we obtain the correct value 

(^) = ^ + ^ + 0(A 2 ). (59) 

Finally, approach II predicts for the average system energy 

1 AO 2 

E s = -d(\n{Z/Z B })/d(3 =— + —f3 + 0(A 2 ). (60) 



It is noted that Eqs. (160]) and (1591) have been derived in [13| using a different method. 
The results obtained for the average system energy via the different approaches, Eqs. (1581) . 
( |59l) . and ( l60i) . are all different, thus providing a nice example of how noncritical use of the 
standard thermodynamic equations (1231) can lead to ambiguous results. 



2. A quantum harmonic dumbbell 



As a final example we consider the harmonic model of a dumbbell coupled to a bath 



quantum mechanically 



36 



371 ] . The thermodynamics of the quantum harmonic dumbbell 



can be constructed within the Wigner distribution method employing the general equations 
( I52al) -( l53l) . The total classical dumbbell+bath Hamiltonian H in Eq. ( I52al) is given by Eqs. 
(I44al) and (T4"5l) . Further, H is inserted into Eq. (1531) to derive the quantum correction 
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The so-obtained expression is, however, quite cumbersome and is not presented. Here, we 
only give the results necessary for the discussion of the thermodynamics quantities. 

The reduced Wigner distribution of the system is given by the general formula Eq. (1541) . 
where H$ is defined via Eq. f)44al) and the bath induced operator reads 

A, = A° 5 - ^A*. (61) 
Here we have introduced the notation 

MO 2 

A = ^A (Xi _ X2)2? (62a) 

a* = M(n 2 + nmx, - x 2 f + ^(n? + nl) + + n|) - ^n 2 , (62b) 

1 N B l N B 

3=1 3=1 

and the frequency Q 2 A of the bath- induced harmonic potential is defined via Eq. (|46|) . The 
bath-induced operator A s is explicitly coordinate-dependent (due to the classical contribu- 
tion (I62ap ) as well as momentum and temperature-dependent (due to the quantum correction 
(]62bp ). The corresponding partition function reads 

Z S = Zg |l + ^ (4fi 2 5 + 2nl + ^ 2 )}, (63) 

where the classical system partition function Z s l is given by Eq. (I47bl) . 

The mean energy of the dumbbell calculated via approach I (employing Eq. fl54"j) ) is given 
by the expression 

(H s ) = (H s ) d + X(H s ) q + 0(X 2 ). (64) 
Here the first term is the classical contribution, 

<*«>- 4 ( i + 1 A) (65) 

and the quantum correction reads 



(66) 



In the limit Q 2 S = 0, Vt\ = Q 2 the dumbbell reduces to two noninteracting Brownian particles. 
Correspondingly, the mean energy Eq. (164")) gives twice of what is predicted by Eq. (|59|) . The 
same is true for the bath induced interaction operators As (cf. Eqs. flST]) and fl55l) ). 
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On the other hand, if we follow approach II, we obtain for the ratio of the partition 
functions of the overall system and the bath 

ZjZ B = Z<§ { 1 - ^ (20 2 s + n? + n*) } (67) 

and correspondingly as prediction for the mean energy of the system 

- d{\n{Z/Z B })/d(3 = A + M ( 2 fl| + n? + ^2) + 0(A2) _ (6g) 

The comparison with Eq. fl64l) shows that approaches I and II give not only different classical 
contributions to the mean system energy, but also very different quantum corrections. It is 
also interesting to note that the strength of the bath-induced potential, does not enter 
the expression for the mean energy, Eq. ([BSD j w hile the mean energy calculated via approach 
I (Eq. fl64|) ) depends sensitively on this quantity. 



V. DISCUSSION AND CONCLUSIONS 

The results presented above demonstrate that the two different approaches to describe the 
thermodynamics of a subsystem can predict very different results if the system-bath coupling 
is not weak. This was already shown earlier for a quantum harmonic oscillator [l^ and for a 



quantum Brownian particle [13j bilinearly coupled to a harmonic bath. The results obtained 
here corroborate and extend these earlier findings. The study also shows that ambiguities 
in the description of reduced thermodynamics already occur in the classical case for more 
complex systems, such as an anharmonic bath (Sec. IIVB II) or if the system under study 
consists of more than a single point particle (Sec. IIVB2p . Based on the results above we 
shall now analyze the two approaches and discuss their advantages and shortcomings. 

The different predictions of the two approaches for the thermodynamic quantities can be 
related to the ambiguities in the definition of the energy of a system that is coupled to the 
environment. If the system-bath coupling is not negligible, the system-bath interaction Hsb 
(or a certain part of it) may be included in the system energy [35( thus resulting in a variety 
of definitions. Furthermore, the results for the thermodynamics quantities may depend on 
the particular physical quantity that is considered fundamental (e.g. the system energy or 
the partition function). 

Approach I associates the mean system energy E$ with the expectation value of the 
system Hamiltonian (Hs)- The definition Eg = (Hs) clearly associates the observable (Eg) 
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with the corresponding physical operator Hs- This definition appears natural and fits into 
the general scheme of the statistical thermodynamics and probability theory. For example, if 
we think of E$ as the mean internal energy (e.g., as of the internal energy of a molecule with 
several vibrational degrees of freedom), the choice E$ = (Hs) is well physically justified. As 
a result of the definition E$ = (Hs), the bath and the system-bath coupling influence the 
mean system energy Eg only indirectly, through the reduced distribution (density matrix) 
of the system ps given by Eq. ffTBl . It is important to note that ps does not coincide with 
the canonical distribution for the isolated system. As a consequence, the thermodynamics of 
the reduced system is described by the relations (|23l)-(l271). Once the definition E$ = (Hs) 
is accepted, no other assumptions are necessary to construct the system thermodynamics. 



Corroborating the results obtained in Refs. 



12 



13 



the present extended study shows that 



the use of this definition gives rise to physically and logically consistent results for both 
quantum mechanical and classical systems. 

The fundamental quantity of approach II, on the other hand, is the partition function 
of the system, Zs- To obtain thermodynamics quantities, approach II involves two steps: 
First, the partition function of the system is identified as the ratio of the total and bath 
partition functions Z/Zb (Eq. ffTBl ). In addition to this choice, approach II assumes that 
the standard thermodynamical relations given by Eqs. ©-(JH]) can be used to calculate 
thermodynamics quantities such as the mean energy, the entropy, and the specific heat. The 
choice Zs = Z/Zb for the partition function of the system appears to be reasonable, notably 
in the limit of weak system-bath coupling. However, as discussed in Sec. IIII CI the partition 
function Zs = Z/Zb corresponds the reduced distribution ps given in Eqs. ( fT3l) . (Ell), which 
does not coincide with the corresponding canonical distribution for the isolated system alone. 
Instead, ps contains an additional bath-induced operator As which, in general, depends on 
the temperature and on the degrees of freedom of the system. In such has been 

shown in Sec. IIII Al the thermodynamic relations for a subensemble, Eqs. (l23l-(!27l). should 
be employed instead of the standard thermodynamic relations (|71)- (|T0l) . provided the system- 
bath coupling is not small. Therefore, there is no any a priori theoretical justification 
for plugging Z s = Z/Z B into Eqs. ff?l)- ffT0]) . and we have to additionally postulate that 
differentiations of F s = — \n(Z/Z B )/ ft give us, according to Eqs. (JTJ)-([9]) , the mean system 
energy, entropy, and the specific heat. Thus the two fundamental assumptions of approach 
II, the choice for Zs and the validity of the standard thermodynamic relations (171)-( fT0l) even 
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for the subensemble, can not be proven within the approach itself. 

It also worthwhile to mention another peculiarity of approach II. As has been pointed 
out in Refs. 
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the mean energy of the system obtained within approach II corresponds 



to the definition 

E s = (H) - (H B ) B = (H s ) + (H B + H SB ) - (H B ) B , (69) 

where (■■■) B denotes averaging over the bath distribution (Eq. (|T5l)). For nonvanishing 
system-bath interaction, the term (H B + Hs B ) — {H B ) B gives an additional contribution to 
E$ that is not present in approach I. As a consequence of the structure of this additional 
term, it is not possible to introduce an operator of the mean energy, whose average will give 
Es- Furthermore, considering for example Es as the mean internal energy of a molecule 
with several vibrational degrees of freedom, it does not appear to be consistent that the 
definition ( 1691) contains contributions which are explicitly determined by the bath degrees 
of freedom. 

To summarize, although approaches I and II give identical result for the thermodynamics 
of a subsystem if the system-bath coupling is negligible, their predictions differ significantly 
for finite system bath coupling. These differences arise because different quantities are 
considered as fundamental in the two approaches and be related to the different definitions 
of the mean energy of the system used. The results obtained above and those presented 
earlier [3] suggest that approach I is superior both from the physical and the logical point 
of view. 
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